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Abstract — In this letter, we propose the resonant tunneling 
diode (RTD) based on a double-barrier graphene/boron nitride 
(BN) heterostructure as device suitable to take advantage of 
the elaboration of atomic sheets containing different domains 
of BN and C phases within a hexagonal lattice. The device 
operation and performance are investigated by means of a self- 
consistent model within the non-equilibrium Green's function 
formalism on a tight-binding Hamiltonian. This RTD exhibits 
a negative differential conductance effect which involves the 
resonant tunneling through both the electron and hole bound 
states of the graphene quantum well. It is shown that the peak- 
to-valley ratio can reach the value of 4 at room temperature for 
gapless graphene and the value of 13 for a bandgap of 50 meV. 

Index Terms — Graphene, Boron nitride substrate, Resonant 
tunneling diode. 



I. Introduction 

THE benefits of the extraordinary intrinsic transport prop- 
erties of graphene [lj are usually hindered by the defects 
of the supporting insulating substrate. Indeed, graphene on 
various substrates such as S1O2 1 2], SiC [3] or other high- 
k insulators El shows a reduced carrier mobility, which is 
usually attributed to surface roughness, charge surface states 
or surface optical phonons. Recently, it has been shown that 
graphene reported on hexagonal boron nitride (h-BN) has 
higher mobility than on any other substrate 0, 0, Q. A 
mobility of 275 000 cm 2 /Vs at low temperature and 125 000 
cm 2 /Vs at room temperature, i.e. as high as for suspended 
graphene, has even been reported 0. It is due to the fact that 
the surface of h-BN is flat, with a low density of charged 
impurities, does not have dangling bonds and is relatively 
inert 0. Hence, with the same atomic structure as graphene 
and a 1.8% higher lattice constant h-BN is becoming a 
very promising candidate as high bandgap insulating material 
to make possible reaching the ballistic transport regime in 
graphene devices even at room temperature and exploiting the 
peculiarities of graphene properties inherent in the massless 
and chiral character of charge carriers. Additionally, it has 
been shown that it is possible to obtain atomic sheets con- 
taining different domains of h-BN and C phases with various 
compositions fTOll , which offers the perspective of designing 
new heterostructures and devices ifTTIl . 
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Fig. 1. Simulated device wherein the graphene sheets are deposited on a 
h-BN substrate and two barriers are thus formed by the BN layers between 
the graphene zones. 

In the field of quantum devices, graphene resonant tun- 
nelling diodes (RTDs) have been proposed and investi- 
gated previously lfT2l . fT3ll . However, while the operation of 
graphene nanoribbon (GNR) RTDs is strongly dependent on 
the device shape lfT2l . the 2D graphene structure studied in 
fT3l assumes that a bandgap opens in graphene epitaxially 
grown on SiC, which is a debated question, in addition to 
the substrate-induced mobility reduction. In this work, we 
thus propose a new graphene RTD based on graphene/h-BN 
heterostructures, where the two barriers are assumed to be 
formed by large bandgap h-BN zones, as schematized in Fig. 
1. Though it may be useful and possible in case of Bernal 
stacking on BN substrate lfT4l . the bandgap in graphene is not 
necessary for this device to operate correctly. 

To investigate the electrical characteristics of this RTD, we 
have developed a self-consistent simulation code based on the 
Green's function formalism lfl~5l to solve the tight-binding 
Hamiltonian suggested in lfT6l . 

II. Physical model and results 

In the structure schematized in Fig. 1, the BN barriers 
are inserted in a monolayer graphene sheet reported on a 
h-BN substrate of thickness 10 nm and dielectric constant 
e r = 3.5 [5]. We assume the transverse width along OY 
direction to be much larger than the typical device length 
along the transport direction, i.e. a few tens of nanometers, 
so that the extension along the transverse direction can be 
considered through Bloch periodic boundary conditions ifTTl . 
The length of the heavily doped regions at both device ends 
and the BN barrier thickness are assumed to be 17 nm and 1.3 
nm, respectively. Our approach is based on the self-consistent 
solution of the 2D Poisson and Schrodinger equations. The 
nearest-neighbor tight-binding Hamiltonian with parameters 
determined in lfT6l is solved using the Green's function method 
fT5l in the ballistic approximation. In this simulation, the null 
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Fig. 2. (a) Self consistent data of potential energy and (b) corresponding 
electron density at different applied biases. The doped concentration in two 
device ends Np = 5 x 10 16 ra -2 . The vertical-dashed lines indicate the 
h-BN regions. 



Neumann conditions are applied at the 2D-domain boundaries 
(dashed line in Fig. 1). In the case of Bernal stacking, 
graphene reported on a h-BN substrate was shown to exhibit 
a finite bandgap, but this bandgap vanishes when graphene 
is misaligned with respect to h-BN lattice lfT4l . Therefore, 
unless otherwise stated, the zero bandgap of graphene was 
assumed in our study. All simulations were performed at room 
temperature. 

In Fig. 2, the self-consistent solutions of potential and 
electron density are plotted for different biases applied to the 
device with well width D = 4.3 nm and doping concentration 
in access regions Njj = 5 x 10 16 m -2 . The low doping in the 
central zone generates a potential barrier which controls the 
threshold voltage of the device. Differently from the conven- 
tional RTDs |[T8lL Fig. 2(b) shows that the local charge density 
is not symmetrical even at zero bias. This feature can be 
understood as a consequence of the device asymmetry between 
the left graphene section coupled to nitrogen atoms and the 
right section bounded to boron atoms of the h-BN lattice. 
Besides, it is shown that the electron density accumulated in 
the quantum well increases when increasing the bias voltage 
from V to 0.4 V and decreases when further raising the bias 
from 0.4 V to 0.6 V. This feature is nothing but an effect of 
resonant tunneling [18 ] as described below. 

We now go to explore the electrical behavior of the device. 
In Fig. 3, we display the I-V characteristics together with 
a diagram showing the LDOS and the corresponding trans- 
mission coefficient at Vuaa = 0.4 V. Indeed, the diagram in 
Fig. 3(a) shows that the large bandgap (« 3.9 eV |[T6li ) in 
the h-BN barriers generates the quantization of carrier states 
in the graphene quantum well, which gives rise to resonant 
tunneling effects. As a consequence, the I-V curves can exhibit 
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Fig. 4. I-V characteristics (a) for different doped concentrations in two 
contacts and (b) for different well widths D. D = 4.3 nm in (a) and Np = 
5 x 10 16 m- 2 in (b). 



a strong negative different conductance (NDC) in both cases 
of gapless and gapped graphene. In principle, the valley 
current in semiconductor RTDs occurs at a finite bias when 
there is no available state in the emitter region for tunneling 
via the confined states in the quantum well. This feature 
is essentially due to the finite bandgap of the emitter. The 
results obtained here are thus explained as follows. Though 
the energy bandgap of graphene may be actually zero for the 
transverse momentum mode corresponding to the Dirac (or 
K) point (i.e., k y = K y = 27r/3v // 3a c ), a finite energy gap 
Eg(k y ) = 2tcc |l — 2cos (a c k y y/3/2)\ still appears for the 
other modes (see in ref. [19]), where tec = 2.5 eV fT6l and 
a c is the carbon-carbon distance. For instance, E g (k y ) ^0.11 
eV for k y = K y + 0.0087T /a c \/3 as seen in the diagram of 
Fig. 3(a). Therefore, since the total current results from the 
contribution of many transverse momentum modes, the valley 
current occurs even in gapless graphene devices as seen in 
Fig. 3(b) when the tunneling corresponding to finite values 
of \k y — K y \ is suppressed at high bias. Moreover, because 
of the zero (or small) bandgap of graphene, it also shows 
that the current valley in this device is much narrower than 
that in conventional RTDs [18]. Besides, when introducing 
a finite bandgap fT2lL ifTH in graphene, both the valley and 
peak currents are reduced while the peak-to- valley current ratio 
(PVR) increases. The latter reaches 4 and 13 for zero bandgap 
and Eg = 50 meV, respectively. 

Especially, from Fig. 3, we find an interesting feature: 
differently from other RTDs [12], [ 18 ] where only the electron 
states in the quantum well (QW) are resonant, in this device, 
both electron and hole states can contribute to the current, 
as a consequence of the chirality of carriers in graphene. This 
point has been also discussed in l20l and was shown to reduce 
the NDC effect in single potential barrier graphene structures. 
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Fig. 3. (a) Diagram showing the LDOS and the transmission coefficient as a function of energy at Vbias — 0-4 V and (b) I-V characteristics of the simulated 
device (for both cases of zero and finite bandgaps). The figures in (a) are plotted for a transversal momentum k y / K y (momentum at the Dirac point). The 
corresponding electron/hole potential profiles are superimposed on the LDOS (solid lines). The doped concentration Np = 5 x 10 16 ra -2 . 



Practically, as we can understand from Fig. 3(a), the confined 
hole state in the QW makes the tunneling from hole states 
in the emitter to electron states in the collector possible. This 
band-to-band (BTB) process contributes to the rapid increase 
of current beyond Vbias- 0-25 V, in addition to the normal 
resonant tunneling of electrons via the electron state in the 
QW. The vanishing of the latter process is responsible for 
the valley regime in the I-V characteristics. When further 
increasing Vfo as , the BTB tunneling through the first electron 
level in the QW occurs, which leads to the rapid re-increase 
of current and a narrow current valley seen in Fig. 3(b). 

To evaluate the role of device parameters on the RTD oper- 
ation, we display in Figs. 4(a) and (b) the I-V characteristics 
obtained for different doping concentrations N& in the access 
regions and for different well widths D, respectively. In fact, 
when increasing Np, the energy spacing between the Fermi 
level (zero-energy point) and the flat potential energy in the 
device ends at zero bias increases. This results in shifting 
the valley region to high bias. Besides, the enhancement of 
BTB resonant tunneling at high bias makes both the peak 
and valley currents higher when increasing Np. It finally 
reduces the PVR which reaches ~ 3.9, 3.8, 2.2 and 1.8 
for N D = 2.5 x 10 16 , 5.0 x 10 16 , 7.5 x 10 16 , and 10 17 
m -2 , respectively. Fig. 4(b) shows that when increasing the 
well width D, the current valley moves to the low bias 
while the NDC effect generally weakens. This feature can 
be straightforwardly understood as the influence of D on the 
position and number of energy levels in the QW. In the case of 
large D (e.g., 6.4 nm in Fig. 4(b)), the multi-level contribution 
gives rise to a quite complex I-V curve with low PVR. 

Finally, though not shown here, the overall current appears 
to be strongly reduced when increasing the BN barrier thick- 
ness, as a consequence of the large bandgap of h-BN. It 
therefore suggests that to observe well the effects discussed 
above, a barrier thickness smaller than 2 nm is mandatory. 

III. Conclusion 

A simulation study of graphene RTDs based on graphene/h- 
BN double-barrier structures has been performed by means of 



the self-consistent solution of the 2D Poisson and Schrdinger 
equations within the tight-binding. It was shown that the 
resonant tunneling and the resulting NDC behavior may in- 
volve both the electron and hole bound states of the graphene 
quantum well. The sensitivity of the electrical characteristics 
to the device parameters as the doping density of access zones, 
the QW thickness and the barrier thickness has been analyzed. 
Though the chiral band-to-band tunneling tends to reduce the 
width of the valley region in the I-V characteristics and to 
induce a rapid re-increase of current, the PVR of NDC effect 
can reach the value of 4 at room temperature for gapless 
graphene and the value of 13 for a small bandgap of 50 
meV which may result from Bernal stacking of graphene on 
h-BN. This work suggests that the engineering of graphene/h- 
BN structures is opening a new route for high-performance 
graphene devices. 
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